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Abstract 

A new transport equation for the intermittency 
factor was proposed to predict separated and 
transitional boundary layers under low-pressure 
turbine airfoil conditions. The intermittent behavior 
of the transitional flows is taken into account and 
incorporated into computations by modifying the 
eddy viscosity, with the intermittency factor, 
7 . Turbulent quantities are predicted by using 
Menter's two-equation turbulence model (SST). The 
intermittency factor is obtained from a transport 
equation model, which not only can reproduce the 
experimentally observed streamwise variation of the 
intermittency in the transition zone, but also can 
provide a realistic cross-stream variation of the 
intermittency profile. 

In this paper, the intermittency model is used 
to predict a recent separated and transitional 
boundary layer experiment under low pressure 
turbine airfoil conditions. The experiment provides 
detailed measurements of velocity, turbulent kinetic 
energy and intermittency profiles for a number 
of Reynolds numbers and freestream turbulent 
intensity conditions and is suitable for validation 
purposes. Detailed comparisons of computational 
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results with experimental data are presented and 
good agreements between the experiments and 
predictions are obtained. 

Nomenclature 

C p pressure coefficient, I — ( U e /U re f ) 2 
FSTI freestream turbulence intensity(%) 

K t flow acceleration parameter, ( u/U 2 )(dU/ds ) 

k turbulent kinetic energy 

L nominal suction surface wetted length 

N non-dimensional spot breakdown rate parameter, 

ncrOl/v 

n spot generation rate 

p static pressure 

Ptotai total pressure 

Re Reynolds number, LU T % i/v 

Re st (fit “ s s )U e /iA 

Reet d t U e /v 

s wetted streamwise distance along suction surface 

Tu turbulence intensity (%), u f fU 

U boundary layer streamwise velocity 

U e local freestream velocity 

U Te f nominal exit freestream velocity 

u T friction velocity 

W magnitude of vorticity 

y n distance normal to the wall 

y + ynU T /v 

7 intermittency factor 

6 momentum thickness 

Xe pressure gradient parameter, (6 2 /u)(dU/ds) 

p molecular viscosity 
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fj, t eddy viscosity 

V /i/p 

Mt Ip 

p density 

a spot propagation parameter 

subscripts 

e freestream 

s onset of separation 

t onset of transition 


1 Introduction 

Flow transition plays an important role in 
turbomachinery applications. Majority of boundary 
layer flows in turbomachines involve flow transition 
under the effects of freestream turbulence, diverse 
pressure gradients, wide range of Reynolds numbers, 
flow separation, and unsteady wake-boundary layer 
interactions. 

Prediction of this type of complex flows is an 
important element in analysis and performance 
evaluation of gas turbine engine components and 
ultimately in the design of more efficient jet engines. 
Especially, in low pressure turbine applications 
prediction of transition behavior is even more 
important for reasons of design efficiency. For low 
pressure turbines the flow is mostly turbulent at 
the high Reynolds number conditions encountered at 
take off and the efficiency is at its design maximum. 
However, due to decrease of Reynolds number caused 
by high altitudes conditions at cruise speeds, design 
based on the sea level conditions tends to under- 
predict -losses and thus leads to substantial drops in 
efficiency (Mayle 1991; Rivir 1996; Lake et al. 2000). 
These losses are attributed to flow separation on 
the suction surface of the turbine blades. At low 
Reynolds numbers with low freestream turbulence, 
the boundary layers on the airfoil surface have a 
tendency to remain laminar and hence the flow may 
separate before it becomes turbulent. This may cause 
a drop in efficiency and result in increase of fuel 
consumption. The impact of such losses is directly 
felt on the operation costs. It has been estimated that 
a 1% improvement in the efficiency of a low pressure 
turbine would result in a saving of $52,000 per year 
on a typical airliner (Wisler, 1998). 

In order to calculate the losses and heat transfer 
on various components of gas turbine engines, and to 
be able to improve component efficiencies and reduce 
losses through better designs, accurate prediction 
of development of transitional boundary layers is 


essential (Mayle, 1991). For an accurate prediction 
of transitional flows under the diverse conditions 
encountered in turbomachinery applications, the aim 
is to And a model providing the following properties: 

• Physically accurate and versatile: the model 
should be capable of accurately predicting 
transitional flows under the diverse conditions 
encountered in turbomachinery applications, 
such as pressure gradients, freestream 
turbulence, wide range of Reynolds numbers, 
unsteady wake-boundary layer interactions and 
flow separation. 

• Computationally efficient and inexpensive: 
The model should not involve extensive 
computational effort and should be relatively 
cheap to compute. 

• Compatible with current CFD methods: The 
model should be easy to implement into existing 
CFD codes without requiring extensive changes 
in computational strategy. 

One of the current methods for predicting flow 
transition is to use the stability theory. In this 
method, stability equations are solved at streamwise 
stations in order to predict the onset of transition. 
This method requires prior solution of the mean flow 
field and returns only the onset point of transition 
without any information on the turbulent part of the 
flow. Another method is using empirical correlations 
in the form of e n . This type of methods also 
require prior solution of the mean flow field. These 
two methods are not compatible with current CFD 
methods. 

One method compatible with current CFD 
methods is the use of low-Reynolds number 
turbulence models. Savill (1993a, 1993b), has 
organized a number of workshops to assess the 
capability of current turbulence models in predicting 
flow transition. The comparisons showed that none 
of the existing models was adequate to predict 
flow transition for a range of flow conditions. 
Westin and Henkes (1997) have also tested a large 
variety of turbulence models and compared models’ 
performances in predicting a few of the T3-series 
transition flow experiments (Savill, 1993a). They 
indicated that none of the models could predict both 
onset location and length of transition for a range of 
flow conditions. This outcome is not at all surprising 
since most of the current turbulence models are not 
designed to predict flow transition. 

An alternative method to this approach is to 
use the concept of intermittency to blend the flow 
from the laminar to the turbulent regions. This 
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approach, although highly empirical, has shown some 
successes in predicting transition behavior. Dhawan 
and Narasimha (1958) correlated the experimental 
data and proposed a generalized intermittency 
distribution function across flow transition. The 
correlation was later improved by Gostelow et 
al. (1994) for flows with pressure gradients subject 
to a range of freestream turbulence intensities. 

Solomon et al. (1995), following the work of Chen 
and Thyson (1971), developed an improved method 
to predict transitional flows involving changes in 
pressure gradients. In this model, the effects 
of changing strearnwise pressure gradient on the 
breakdown physics and spot spreading rates are 
taken into account. This is accomplished by varying 
the spot spreading angle and propagation parameter 
through the transition zone according to the local 
pressure gradient parameter. 

Steelant and Dick (1996) proposed a transport 
equation for intermittency, in which the source 
term of the equation is developed such that the 7 
distribution of Dhawan and Narasimha (1958) across 
the transition region can be reproduced. Steelant 
and Dick used their model, coupled with two sets 
of conditioned Navier-Stokes equations, to predict 
transitional flows with zero, favorable, and adverse 
pressure gradients. However, since their technique 
involved the solution of two sets of strongly coupled 
equations, the method is not compatible with existing 
CFD codes, in which only one set of Navier-Stokes 
equations is involved. Moreover, the model was 
designed to provide a realistic strearnwise 7 behavior 
but with no consideration of the variation of 7 in the 
cross-stream direction. 

Cho and Chung (1992) developed a k — 
e - 7 turbulence model for free shear flows. 
Their turbulence model explicitly incorporates the 
intermittency effect into the conventional k — e model 
equations by introducing an additional transport 
equation for 7 . They applied this model to compute 
a plane jet, a round jet, a plane far wake, and a plane 
mixing layer with good agreements. Although this 
method was not designed to reproduce flow transition 
it provided a realistic profile of 7 in the cross-stream 
direction. 

Suzen and Huang (1999) improved the 
intermittency transport equation by combining 
the best properties of Steelant and Dick’s model and 
Cho and Chung’s model. Their model reproduces 
the strearnwise intermittency distribution of Dhawan 
and Narasimha (1958) and also produces a realistic 
variation of intermittency in the cross-stream 
direction. The model is capable of predicting flow 
transition under diverse conditions. The predicting 


capabilities of this model have bevi validated against 
T3- series experiments of Seuill (1993a, 1993b) 
and low-pressure turbine experiments of Simon et 
al. (2000) with good success (Suzen and Huang. 
1999, 2000; Suzen et al. 2000). 

In the current research we focus on a further 
validation of the intermittency transport model of 
Suzen and Huang (1999) against the experiments 
of Hultgren and Volino (2000J Hultgren and 
Volino’s experiments investigated the effects of 
freestream turbulence and Reynolds number on 
separated and transitional boundary layers under 
low pressure turbine airfoil conditions. In their 
experiments, a flat plate boundary layer subject 
to a strearnwise pressure gradient was studied. 
The superimposed pressure gradients were produced 
by attaching a two-dimensional contoured shape 
to the wall opposite to the tesc surface and by 
applying suction on the contoured wall. The 
resultant pressure profile represents that on the 
suction side of the Pak-B airfoi . The experiments 
covered a range of flow conditions including Reynolds 
numbers between 50,000 and 300 TOO and freestream 
turbulence intensities between 0 and 7%. These 
cases cover a realistic range of operating conditions 
from takeoff to cruise. On the tet wall, velocity and 
turbulent kinetic energy and irrvrmittency profiles 
were measured at fourteen streemwise stations. In 
the measurements, quantities si c .1 as skin friction 
coefficients, transition start and end locations, and 
the locations of separation and 1 eattachment were 
also determined. Further details of the measurements 
and experimental data were given by Volino and 
Hultgren (2000). 

The experiments of Hultgren and Volino (2000) 
provide a good set of data fo the development 
and validation of the models tor flow transition. 
In the next section, the intermittency transport 
model is presented and implementation of the model 
is described along with the empirical correlations 
employed for the onset of transitior. In Section 3, the 
numerical details of the prediction orocess are given. 
In Section 4, the predictions of tin new intermittency 
model are compared against the e cperimental data. 
Finally, conclusions are provided n Section 5. 

2 Transport Tv odd 
for the Interm i t tencv 

In this section, the transport model for 
intermittency is presented. The model combines 
the transport equation models of Steelant and Dick 
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(1996) and Cho and Chung (1992). Details of the 
development and implementation of the transport 
model are given :n Suzen and Huang (1999,2000), 
Suzen et al. (2000). 

The model equ ation is given by: 


dpi dpu 3 7 
dt dx : 


(1 - 7) 


(1 - F)2Co Py /w^f(s)f(s) 


+ 

+ 

+ 


pf Cil k 3/ 2 m duj d'y \ 

V k X3 hxj t (ukUk) 1 / 2 dxj dxj ) 

k 2 d'y dj 
6 dxj dxj 

^—(((1 -l)l°ytP+ (1 -l)<7i t pt)^)- ( 1 ) 


The distributed breakdown function, f(s) has the 
form: 


. as' 4 -f 6s' 3 4- cs 2 + ds f + e 

/(s) = ^Th • (2) 

where s' = s - s is the distance along the 
streamline coorcinate, and $ t is the transition 
location. The coefficients are; 

a = b =-0.4906 c = 0.204(^)“° 5 

d = 0.0 e = 0.04444(2 j) -1 ' 5 
h = iOe and g = 50. (3) 


The shear stresses are defined as: 


C 2 — 0.16 and Cz — 0.15. 

The int.ermittency is incorporated into the 
computations simply by multiplying the eddy 
viscosity obtained from a turbulence model, fi t , by 
the inter mitt ency factor, 7. Simon and Stephens 
(1991) showed that by combining the two sets 
of conditioned Navier-Stokes equations and making 
the assumption that the Reynolds stresses in the 
nonturbulent part are negligible, the intermittency 
can be incorporated into the computations by 
using the eddy viscosity, which is obtained by 
multiplying the eddy viscosity from a turbulence 
model, /xt, with the intermittency factor, 7. That 
is, 

Pt = IPt (6) 

is used in the mean flow equations. It must be noted 
that 7 does not appear in the generation term of the 
turbulent kinetic energy equations. 

The intermittency model had been implemented 
into TURCOM code developed by Huang and 
Coakley (1992) and validated against low-pressure 
turbine experiments of Simon et al. (2000) (Suzen 
et al., 2000). 

It is essential that the turbulence model selected 
to obtain pit must produce fully turbulent features 
before transition location in order to allow the 
intermittency to have full control of the transitional 
behavior. Menter’s (1994) SST model satisfies this 
requirement. It produces almost fully turbulent 
flow in the leading edge of the boundary layer and 
therefore it is used as a baseline model to compute pit 
and other turbulent quantities in the computations 
(Suzen et al., 2000). 

The value of no used in evaluating the constants 
given by (3) is provided by the following correlation 
for zero pressure gradient flows (Suzen et al., 2000); 


dui duj 2 duk , 

dxj + dxi 3 dxk %3 



(4) 


The blending function F is constructed using a 
nondimensional parameter, kjWv, where k is the 
turbulent kinetic energy and W is the magnitude of 
the vorticity. The blending function has the form: 


her = (nz/ 2 /t/ 3 )cr = 1.8 x 10 11 Tu^ A . (7) 


When flows are subject to pressure gradients, the 
following correlation is used: 


no _ f J^(l~exp(0.75xl0 6 K t Tu 0 7 )) ^ < Q ^ 

" l 10-3227K?— J Kt > o (8) 


F — tanh 4 


k/Wu 

200(1 — 7°' 1 )°-3 ' 


(5) 


with M defined as: 

M = (850Ttt“ 3 - 100Tu“ 0 - 5 + 120) 


The model constants used in Equation 1 are: 
o lt — o lt — 1.0 Co ~ TO Ci — 1.6 


where, (tio)zpgJs the value for flow at zero pressure 
gradient and can be obtained from equation (7) 
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and K t = (. v/Ut)(dU/dx) t is the flow acceleration 
parameter. The favorable pressure gradient part of 
the above correlation (for K t > 0) is from Steelant 
and Dick (1996). The portion of the correlation for 
adverse pressure gradient flows, K t < 0, is formulated 
using the transition data of Gostelow et al. (1994) and 
Simon et al. (2000) (Suzen et al., 2000). 

The current approach uses the intermittency 
transport model to : obtain the intermittency 
distribution for the transitional flows, while the 
onset of transition is defined by the correlations 
discussed in the next section. 

Attached-Flow Transition 


Abu-Ghannam and Shaw (1980) suggested that 
the onset of transition for attached flows can 
be obtained by correlating the boundary layer 
momentum thickness Reynolds number to the 
freestream turbulence intensity according to: 


Ree t — 163 4- exp 


F(X e ) - 


F{Xe) 

6.91 


Tu 


(9) 


where 


F(\q) = 6.91 -f- 12.75Ag -f 63.64A0 2 for \e < 0, or 

F{\ e ) = 6.91 + 2.48A# - 12.27A*? 2 for A* > 0. 


Although Abu-Ghannam and Shaw correlation 
shows good agreement with experimental data for 
flows with zero and adverse pressure gradients, the 
model is not very sensitive to flows subject to strong 
favorable pressure gradients, in which one would 
expect the transition to be delayed as a result of flow 
acceleration (Suzen et al., 2000). 

To allow for a more sensitive response to strong 
favorable pressure gradients while maintaining the 
good features of Abu-Ghannam and Shaw in adverse 
pressure gradient region, the transition criterion was 
re-correlated to the freestream turbulence intensity, 
Tu, and the acceleration parameter, K t , 

Reet = (120 + 150TiT 2/3 )coth[4(0.3 - K t x 10 5 )], (10) 


where K t was chosen as the maximum absolute value 
of that parameter in the downstream deceleration 
region (Suzen et al., 2000). Equation 10 was designed 
to have a better fit of the available experimental 
data: while the correlation fits the transition data 
well for flows under adverse pressure gradients, it 
was purposely designed to rise rapidly as K t becomes 
positive. This measure is to reflect the fact that 


the flow becomes less likely to have transition when 
subject to favorable pressure gradients. 

Separated-Flow Transition 

Roberts (1980) proposed a semi-empirical theory 
to predict onset of transition within a laminar 
separation bubble over the airfoil suction surface. 
The transition Reynolds number, Re st , which is based 
on the length defined between the onset location of 
separation and that of transition, is correlated as 
a function of a turbulent factor representing effects 
of the external turbulence level and its disturbance 
spectrum. The model was simplified by Davis et 
al. (1987) to only a function of turbulence intensity, 

Re st = 2.5 x 1 0 4 log 10 coth(0. 1732 Tu), (11) 


where Tu is the freestream turbulence intensity value 
at the onset of separation. Although this model was 
originally proposed for swept wing flows, it had been 
used widely for predicting onset of transition in a 
variety of separated flows. However, comparisons 
of recently available data for onset of transition in 
separated flows suggest that the correlation for the 
onset of transition is better represented by a function 
of more than one parameter. 

In Figure 1, Davis et al. (1987) correlation is 
given along with the data of Roberts (1980) (the 
original data used to develop Davis’ correlation), 
Simon et al. (2000), Hultgren and Volino (2000) and 
Lake (1999). Each point (except for the data of 
Robert (1980)) is numbered to indicate the Reynolds 
numbers of the cases based on exit velocity and 
suction surface length. There is some scatter in the 
data especially for high Reynolds number cases of 
Lake (1999). In order to capture the scatter of the 
data, we incorporated the effect of Reynolds number 
into the correlation for Re s t- The new correlation 
expresses Re s t in terms of the turbulence intensity 
(T u ) and the momentum thickness Reynolds number 
at the point of separation (ReQ S ) and is given as; 

Re s t = 874Re$ 5 n exp[— 0.4Tu] (12) 

In Figure 2 the new correlation is compared with 
the experimental data for different values of Re$ s . In 
contrast to Davis et al. (1987) correlation, the current 
formula seems to provide a better representation of 
the data scattering. In the present paper Equation 12 
is used as a replacement of Davis et al. correlation to 
predict the onset of separated-flow transition. 
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3 Numerical Details 

As mentioned earlier, the experiments made 
use of a flow suction on the upper (contoured) 

. surface to arrive at the desired pressure distribution. 
Although the shape of the upper contoured 
surface is known, no information is available for 
the conditions of the suction on this surface. 
This somewhat makes the computation of the 
problem ill-posed. In order to provide the upper- 
boundary of the computational domain, we defined 
a prescribed streamline, obtainable by integrating 
the experimental measured velocity profiles, as the 
upper boundary. Once the streamline is defined as 
the upper computational boundary, a slip boundary 
condition can be applied at the upper surface. 

For example, in Figure 3 the streamlines obtained 
by using different values of 

F = / ^d y, (13) 

Jo 

are presented for the case with Re — 300, 000 and 
FSTI— 7%. In this figure F 0 is the F obtained from 
Equation 13 at the first station taking h as the 
maximum height of the measured y position. To 
define the streamlines, we performed the integration 
of measured velocity profiles based on Equation 13 in 
all x stations such that F is a constant. Any one of 
the streamlines shown in Figure 3 can be used as the 
outer boundary of the computational domain, as long 
as they are well outside the boundary layer. Similar 
streamline curves can be generated for other cases. 

In the computations, we selected a streamline 
corresponding to F = 0.7Fo as the upper boundary 
of the computational domain. The shapes of the 
upper streamlines are presented in Figure 4 for all 
of the cases considered. In this figure the symbols 
show the points obtained from integration of the 
experimental velocity profiles and the lines indicate 
the boundary shapes used in the computations. 
By varying different streamline positions, we have 
found the streamline corresponding to F = 0.7Fo 
is sufficiently remote from the boundary layer to 
be considered as freestream. Non-uniform meshes 
consisting of 140 x 60 grid points are used for all 
calculations. All grids have first y + values less than 
0.5. A grid refinement study for the case with Re — 
50,000 and FSTI=7% (using a 250 x 100 mesh) has 
revealed that the solution obtained by the current 
mesh can be accepted as a grid independent solution. 

Although the last experimental station is located at 
x/L = 1.06, the computational domain is extended 
to x/L — 1.4 in the streamwise direction to . allow 
the outflow boundary condition to be applied in 


the computations. We have extended the upper 
streamline boundary up to x/L — 1.25 and a 
constant pressure outflow boundary conditions is 
applied downstream of that station. Although this 
extension of the domain may appear to be arbitrary, 
it is necessary to avoid having reversed flow close 
to the exit boundary. Since no reverse flow was 
encountered for all cases after x/L ^ 1.25, the 
extension of the domain up to x/L = 1.4 seems 
reasonable. The plate surface is assumed to be 
an adiabatic wall and a constant-pressure outflow 
boundary condition is applied on the exit plane for 
all cases. 

In order to obtain accurate inlet profiles for 
the computations, a laminar computation over a 
fiat plate is performed. From this computation, 
the profile matching the momentum thickness the 
first experimental station is extracted as the inlet 
condition for the computations. 

4 Results and Discussion 

In this section the comparison of the predictions 
with experimental data of Hultgren and Volino (2000) 
are presented. 

The experiment considered in this study involves 
two different Reynolds numbers, Re = 300,000 
and Re = 50, 000 (based on nominal exit velocity, 
U ref : and wetted plate length, L = 0.208m) and 
two freestream turbulence intensities, FSTI— 7% and 
FSTI— 0.2%. Figure 5 depicts the dynamic interplay 
between transition and separation when subject 
to variations in Reynolds number and freestream 
turbulence effects. The first case considered is with 
Re — 300,000 and FSTI=7%. For this case, the 
transition position is at x/L — 0.66 and only a tiny 
flow separation is observed near x/L « 0.7. As the 
Reynolds number is decreased from 300,000 to 50,000 
while keeping the freestream turbulence intensity the 
same, the onset of transition is delayed until x/L = 
0.85. As a result of the delay of the transition, a 
sizable flow separation is encountered between x/L as 
0.7 and x/L ss 1.0. On the other hand, when the 
turbulence intensity is decreased from 7% to 0.2% 
while the Reynolds number is maintained at 300,000, 
the transition onset is delayed until x/L « 0.79. 
As a result, a small but visible separation bubble is 
observed in the vicinity of x/L « 0.75. The last 
case involved a simultaneous decrease of Reynolds 
number and freestream turbulence intensity, from 
300,000 to 50,000 and from 7% to 0.2%, respectively. 
In this case the onset of transition is delayed to 
x/L « 1.03 and a massive separation extending 
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from x/L « 0.7 to x/L as 1.25 is observed. 
The current prediction is an attempt to mimic the 
above-mentioned dynamic behavior of the interaction 
between transition and separation when subject 
to changes in Reynolds number and freestream 
turbulence intensity conditions. 

For the case involving Re = 300, 000 and 

FSTI=7%, experimental onset point of transition 
is at Ree — 336 which corresponds to a location 
x/L — 0.66. However, the onset point of 

transition for this case seems to be lower than other 
experimental data trend compiled by Mayle (1991) 
and Savill (1993a) and also than the value given by 
the correlation of Abu-Ghannam and Shaw (1980). 
For example, the freestream turbulence intensity 
at the point of transition is Tu — 1.7% and the 
acceleration parameter is, K t = —2.5 x 10" 6 , by 
applying (Equation 10), the predicted onset point 
yields a much earlier location (at Ree — 230, 
corresponding to x/L = 0.55). To show the effect 
of the transition positions, we have performed two 
different computations for this case: one using the 
experimental onset point of transition and the other 
using the onset point obtained from the correlation, 
Equation 10. 

The predicted freestream velocity and pressure 
coefficient distributions along the surface are 
compared with experimental data, as shown in 
Figures 6a and 6b, respectively. In these figures, the 
results of the computation using the experimental 
onset of transition point is titled as “Computation 
1” and the one obtained by utilizing Equation 10 is 
denoted as “Computation 2” . Roth results show very 
good agreements with experimental data and only 
slight differences in the profiles shapes near x/L « 
0.65 are observed between the two computations. 

The velocity profiles for the two computations 
are compared with experimental data in Figures 6c 
and 6d. In the experiment, the flow remains 
laminar before station x/L — 0.69. At station 
x/L = 0.69 the velocity profile indicates that 
the flow is on the verge of separation. At the 
next station {x/L = 0.75), the boundary layer 
is attached and the flow continues to develop as 
an attached turbulent boundary layer. Although 
it is not seen from the measured velocity profiles, 
it was reported experimentally that a very small 
separation region exists between the two measured 
stations, x/L — 0.69 and x/L — 0.75. In contrast, 
no separation was observed in both computations 
although the comparison of the velocity profiles shows 
the predictions agree very well with experimental 
data. Both computations produce similar results 
except in the region x/L & 0.7. Due to the fact 


that the onset point of transition of computation 2 is 
slightly upstream, the result of computation 2 is fuller 
and is less likely to separate than that of computation 
1. Otherwise, the differences between the two profiles 
are small. 

The intermittency profiles of the two computations 
are compared with experimental data in Figures 7a 
and 7b. As can be seen from Figure 7a, the flow 
is laminar at the first seven stations. In both the 
experiment and the computation 1, transition begins 
at station x/L — 0.66 and in computation 2 the 
transition starts at x/L — 0.55. While the transition 
of the computation 2 is ahead of the computation 
1, the length of transition is somewhat longer due 
to the fact that a small value of spot generation 
rate, ncr, is calculated. The maximum value of 
the intermittency in the profiles of computation 1 
reaches unity near station x/L = 0.81 whereas 
in computation 2, the intermittency reaches unity 
only after x/L = 0.88. Overall, computation 1 
shows slightly better agreement with experiments 
than computation 2 even though the differences of 
the two results are mainly limited to the near wall 
region. 

In Figures 8a and 8b, comparisons of the turbulent 
kinetic energy profiles are made. As can be seen 
from these figures, the freestream decay of turbulent 
kinetic energy throughout the plate are made to 
match the decay of freestream turbulence. This 
match provides the initial conditions for dissipation 
of the turbulent kinetic energy (see Suzen and 
Huang, 1999, for details). The predicted profiles 
of the turbulent kinetic energy seem to follow the 
experimental trend well even though the values are 
slightly larger before x/L = 0.7 and smaller after 
that. Again, the differences of the two computations 
were only limited to region near x/L = 0.7. 

The next case considered involves the same 
freestream turbulence intensity (FSTI=7%) while 
Reynolds number is reduced from 300,000 to 50,000. 
In the experiment, the flow separated at x/L — 
0.63 and the onset of transition was observed at 
x/L = 0.85. In the computations, separation 
took place at x/L = 0.66 and the onset point of 
transition calculated from Equation 12 was at x/L = 
0.88. The freestream velocity and pressure coefficient 
distributions are compared with experimental data, 
as shown in Figures 9a and 9b, respectively. The 
agreements with the experimental data is good except 
in the region after x/L — 0.95, where it can be seen 
that the freestream velocity distribution is slightly 
over-predicted while the pressure coefficient is slightly 
under-predicted. 

The comparisons of the velocity profiles are given 
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in Figures 9c and 9d. The profiles are in good 
agreement with the experiment up to x/L ~ 0.88. 
The experimental data shows that the flow attaches 
after x/L = 0.94 whereas the predictions indicate 
that the separation is extended up to x/L — 1.0. 

The intermittency profiles are compared to 
experimental data in Figure 10a and 10b. The general 
profile trend is well predicted even though the exact 
profile shapes are not well captured. In Figure 11a 
and lib, comparisons of turbulence kinetic profiles 
are shown. It can be seen that the comparisons 
between the predictions and experimental data are 
favorable. 

The next case involves FSTI=0.2% and Re = 
300,000. In the experiment, the flow separated 
at x/L — 0.67, and the transition onset was 

observed between stations x/L ~ 0.75 and x/L = 
0.81. The prediction shows that separation was at 
x/L — 0.66 and the onset of transition (obtained 
from Equation 12) was at x/L ~ 0.79 . The 
comparisons of the freestream velocity and pressure 
coefficient distributions are presented in Figure 12a 
and Figure 12b, respectively. Both profiles show very 
good agreements with the experimental data. 

The comparison of the velocity profiles are given 
in Figures 12c and 12d. The agreement between 
the predictions and experimental data is very 
good. It should be noted that there seems to be 
some discrepancies between the prediction and the 
measurement in the near-wall region near the flow 
separation region ( x/L ^ 0.8), this difference is 
caused by the failure of the hot wire measurement 
for the flow reversal. 

The intermittency profiles are compared to 
experimental data in Figure 13a and 13b. The trend 
for the streamwise development of the intermittency 
factors seems to be well predicted, even though 
the model predicts a less diffusive behavior of the 
intermittency profiles near the freestream region. 

The turbulence intensity profiles along the surface 
are compared to experimental data in Figure 14a 
and 14b. Although a general trend of the turbulent 
kinetic energy profiles is captured, it can be seen 
that the prediction gives rise to larger values of the 
turbulence kinetic energy for flow up to x/L ^ 0.85. 

Finally, both Reynolds number and the freestream 
turbulence intensity are decreased to 50,000 and to 
0.2%, respectively. For this case, flow separation 
was observed at x/L = 0.63 and the onset of 
transition was observed between stations x/L = 1.0 
and x/L — 1.06 in the experiment. The flow did 
not reattach at the last measured station (x/L=1.06). 
The predictions indicates that the flow separation 
is at a;/L = 0.64 and the onset of transition 


calculated from Equation 12 is at x/L — 1.03. The 
predicted freestream velocity distribution and the 
pressure coefficient distribution are compared with 
experimental data, shown in Figure 15a and 15b, 
respectively. As can be seen the comparisons, the 
agreements between the predictions and experimental 
data are very good. 

The comparisons of the velocity profiles are shown 
in Figures 15c and 15d. As can be seen from 
the figures, good agreement is also observed for 
the velocity profiles. Once again, the discrepancy 
between the prediction and experiments near the wall 
is caused by the failure of the hot wire measurement 
in the separation region. 

The intermittency profiles are compared to 
experimental data in Figure 16a and 16b. As can 
be seen from the figures, the flow remains laminar up 
to station x/L = 1.0. At the last measure station, 
x/L — 1.06, the predicted intermittency factors show 
a larger magnitude than the experiments. 

The turbulence intensity profiles are compared to 
experimental data in Figure 17a and 17b. It can be 
seen that the prediction gives rise to larger values 
of the turbulent kinetic energy profiles for region 
x/L > 0.65. Even though the agreement is not that 
good, due to the fact that the transition only occurs 
at x/L = 1.03 its impact to the velocity profiles is 
not that pronounced. 

5 Concluding Remarks 

A new transport equation for the intermittency 
factor is employed to predict a recent transitional 
boundary layer flow experiment under low pressure 
turbine airfoil conditions. The intermittent behavior 
of the transitional flows is taken into account by 
modifying the eddy viscosity with the intermittency 
factor. The new transport model not only can 
reproduce the experimentally observed streamwise 
variation of the intermittency in the transition zone, 
it also provides realistic cross-stream variation of the 
intermittency profile. Computations are performed 
for two different Reynolds numbers and two different 
values of free stream turbulent intensities. Detailed 
comparisons with experiments are made for pressure 
coefficients, velocity, intermittency and turbulent 
kinetic energy profiles. Overall, good agreement 
with the experimental data is obtained. It has 
been demonstrated that the predictions accurately 
mimic the dynamic behavior of the interplay 
between the transition and separation when subject 
to variations of Reynolds number and freestream 
turbulent intensity conditions. 
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Figure 1: Onset of separated flow transition. 



Tl (%) 


Figure 2: Onset of separated flow transition with new correlation (Equation 12). 
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Figure 5: Overall view of effects of Reynolds number and freestream turbulence intensity. 
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7: Comparison of intermittency profiles along the surface- Re = 300, 000 FSTI=7%. 




Figure 8: Comparison of turbulence intensity profiles along the surface- Re = 300, 000, FSTI=7%. 
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Figure 9: Prediction of experiment of Hultgren and Volino (2000)- Re — 50,000, FSTI=7% 
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Figure 10: Comparison of intermittency profiles along the surface- Re = 50, 000, FSTI— 7%. 




Figure 11: Comparison of turbulence intensity profiles along the surface- Re — 50,000, FSTI=7%. 


American Institute of Aeronautics and Astronautics 
















c)2001 American Institute of Aeronautics & Astronautics or Published with Permission of Author(s) and/or Author(s)’ Sponsoring Organization 



20 

American Institute of Aeronautics and Astronautics 









c)2001 American Institute of Aeronautics & Astronautics or Published with Permission of Author(s) and/or Author(s)' Sponsoring Organization. 



3 C 

x/L=0.28 0.33 

0.39 

0.45 

0.51, 

0.57 

0.63 

j 

i 

> <i 

' i 

• < 


c 

■ experiment 
) calculation < 



(a) 


0 1.2 2.4 3.6 4.8 6 7.2 8.4 

T 



Figure 16: Comparison of interinittency profiles along the surface- Re — 50,000, FSTI— 0.2%. 
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Figure 17: Comparison of turbulence intensity profiles along the surface- Re = 50,000, FSTI=0.2%. 
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